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Introduction 


Wave  propagation  along  a perfectly  conducting  cylindrical  wire 

parallel  to  a plane  interface  between  two  dissimilar  materials  has  been 

an  area  of  theoretical  and  practical  interest  for  several  decades.  During 

this  time  the  problem  has  been  treated  in  varying  degrees  of  approximation 

by  many  investigators.  Among  the  more  relevant  to  this  work  are  the  treatments 
1 2 

of  Wait  and  Chang  and  Olsen  . Wait  postulated  a primary  longitudinal  Hertz 
potential  which  generated  a cylindrically  symmetric  TM  field  surrounding  the 
wire.  He  then  expanded  the  field  into  a spectrum  of  plane  waves  and  scattered 
these  plane  waves  from  the  interface  to  obtain  the  total  electric  field.  Requir- 
ing the  longitudinal  component  of  this  total  electric  field  to  be  zero  on 
the  surface  of  the  wire  he  obtained  an  equation  for  the  modal  propagation 
constants.  Chang  and  Olsen  studied  this  modal  equation  in  detail  and  showed 
the  existence  of  not  only  the  usual  transmission  line  mode  but  also  a so- 
called  "earth-attached"  or  "fast- wave"  mode.  However,  all  of  this  work  was 
done  in  the  thin-wire  approximation.  That  is,  the  only  current  accounted 
for  on  the  wire  was  the  angular  average  longitudinal  electric  current. 

If  the  wire  is  not  thin  and/or  is  less  than  a fraction  of  a wavelength 
away  from  the  interface,  one  would  expect  that  the  longitudinal  current  would 
have  some  azimuthal  variation  and  that  there  would  be  some  azimuthally 
directed  currents.  Thus,  one  is  led  to  inquire  as  to  the  conditions  under 
which  these  angularly  dependent  currents  can  be  neglected  in  formulating 
the  modal  equation. 

In  addressing  ourselves  to  the  question  posed  above,  we  formulate  the 
problem  in  general  and  then  express  the  modal  equation  retaining  only  the 
zero-order  and  first-order  terms  in  the  azimuthal  variation.  Then  by 
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comparing  the  first- order  terms  with  the  zero-order  ones,  we  obtain  a 

criterion  for  determining  the  error  incurred  in  retaining  only  the  zero 

order  terms.  Some  work  in  this  regard  was  carried  out  by  Grinberg  and 
3 

Bonshtedt  but  their  results  appear  to  be  limited  to  cases  where  the 
propagation  constant  is  unmodified  by  the  presence  of  the  first  order  currents. 
That  is,  they  have  assumed  conditions  such  that  the  zero-order  equations 
de-couple  from  the  higher  order  ones.  The  approach  presented  here  takes  full 
account  of  the  coupling. 


Formulation 

The  geometry  to  be  analyzed  is  shown  in  Figure  1.  The  wire  is  infinite 

ik.az  • t 

in  the  z-direction  and  all  field  quantities  are  assumed  to  vary  as  e e"  w 
where  kj  = • Similarly  k2  = w/vjTe^"  and  may  be  complex  to  account 

for  loss  in  medium  2.  A cylindrical  coordinate  system  is  shown  with  its  axis 
coincident  with  the  axis  of  the  wire.  We  postulate  an  infinite  sum  of  two 
dimensional  electric  and  magnetic  multipoles  of  unknown  amplitudes  located 
on  the  axis  of  the  wire.  Expressing  the  field  of  each  multipole  as  a spectrum 
of  plane  waves  and  scattering  each  plane  wave  from  the  interface  one  could 
obtain  an  expression  for  the  total  tangential  electric  field  exterior  to 
the  wire.  Then,  requiring  the  tangential  field  components  to  be  zero  on  the 
wire  surface  one  would  obtain  an  infinite  set  of  simultaneous  linear  equations 
for  the  amplitude  of  the  multipoles.  Being  homogeneous,  these  equations 
require  for  non-trivial  solution  that  the  determinant  of  the  coefficients  be 
zero.  This  requirement  would  yield  the  modal  equation. 

Before  proceeding  with  this  formulation,  however,  we  first  divide  the 
modes  into  three  categories  of  which  only  one  will  be  studied.  The  categories 
are:  those  with  no  zero-order  magnetic  multipole  field,  those  with  no  zero 
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order  (angularly  independent)  electric  multipole  field,  and  hybrid  modes 
with  both  zero-order  multipole  fields.  In  this  treatment  we  will  concern 
ourselves  only  with  the  modes  of  the  first  category  because  among  them  are 
the  modes  studied  in  the  thin-wire  approximation  by  previous  investigators. 
Moreover,  we  will  treat  only  modes  which  are  symmetric  about  the  x-z  plane. 
Thus,  we  postulate  the  presence  of  a zero-order  electric  multipole  and  first 
order  electric  and  magnetic  multipoles  and  assume  that  all  higher  order 
multipoles  have  negligible  amplitude. 

The  postulated  multipoles  will  be  represented  by  their  Hertz  vector 
potentials  as  follows.  The  zero-order  electric  multipole  potential  is 


n 

o 


K1 


(1) 


where 


= (1  - ex2)* 


and  Ixn  CC)  > 0 


Now,  we  note  that  one  may  generate  from  this  potential  the  potential  for  the 
first-order  electric  multipole  by  merely  differentiating  with  respect  to  the 
distance  d from  the  multipole  to  the  interface  and  dividing  by 
That  is, 

Q 
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If  there  were  to  be  a zero-order  magnetic  multipole,  its  potential  would  be 


(CkjP) 


(3) 


Now,  the  symmetric  first  order  magnetic  multipole  field  corresponds  to  the 
antisymmetric  first  order  potential.  This  potential  can  be  obtained  from 
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(3)  by  taking  the  negative  of  its  derivative  with  respect  to  y and 


dividing  by  kjC  . Thus 


’I  - h 'o ' jr  wV,si"  ♦ (4) 

Next,  each  of  the  potentials  (1),  (2),  and  (4)  will  be  used  as  a primary 
source  a distance,  d,  from  the  interface  and  in  each  case  the  components  of 
the  total  exterior  electric  fields  tangent  to  the  wire  surface  will  be 


derived. 


Considering  first  the  zero-order  electric  multipole  field,  we  write 
its  potential  as  a plane  wave  spectrum  in  the  form 


^e  1.1  f 1 
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, -u  |k  x-D|  -iXk  y 
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where  u^  = [X  -x,  ] and  D = k^d.  Now,  following  a procedure  similar  to 
that  outlined  by  Wait1  we  introduce  transmitted  and  reflected  Hertzian 
potentials  of  electric  and  magnetic  types  as  follows. 


, 1 *B(M  -iXk  y 
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for  x > 0 (6) 
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for  x < 0 (8) 


for  x < 0 (9) 


where  u_  = [X  +a  -n  1 , n = 2/k,  and  n and  it.,  are  magnetic  Hertz  vectors 
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Now,  computing  the  tangential  components  of  electric  and  magnetic  fields  from 


these  potentials  and  requiring  continuity  at  x=0  we  obtain, 
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1+R(x)  = ^2 
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Thus,  the  z component  of  the  total  electric  field  for  x > 0 is. 
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where  p'  = [ (x+d)2+y2] *^2  , and  D = k^d. 


Turning  next  to  the  first  order  electric  potential  (2) , recalling  that 
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-*e  1 ""o,  and  noting  the  linearity  of  the  reflection  piocess  we  see 
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where  p'  and  <J>'  are  cylindrical  coordinates  referred  to  x = -d,  y = 0;  i.e., 

x+d  x+d 


cos<J> ' = 


p'  r.  ,,.2  2,1/2 

[(x+d)  +y  ] 


(15) 


and 


sin<t>'  = — Y-  = 


p'  [ (x+d)2+y2] 1/2 


(16) 
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Similarly,  by  first  computing  tt^,  ir^R,  and  and  then  substituting  in. 
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we  obtain. 
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In  order  to  obtain  fields  due  to  a first-order  magnetic  multipole,  we 
begin  with  a zero-order  magnetic  multipole  in  anticipation  of  obtaining  the 
desired  result  by  differentiation.  Thus  we  proceed  to  express  the  zero-order 
potential  in  terms  of  its  plane  wave  spectrum  as, 


2 
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As  before  we  postulate  the  following  scattered  potentials. 
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where  we  emphasize  that  “V^  and  ~tt^n  are  electric  Hertz  potentials.  Imposing 


continuity  of  tangential  fields  at  the  interface  yields. 
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1 + R*(X)  = 1 ■ — u. 
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Thus,  the  z component  of  the  total  electric  field  for  x > 0 is  generated  by 
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jtqM  . In  particular. 
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and  by  differentiating  with  respect  to  y and  dividing  by  kj  we  obtain, 
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as  the  z- component  of  the  total  electric  field  for  x > 0 due  to  a first  order 
magnetic  multipole  source.  Likewise, 


_ m 3 m . 3 .m  m. 

E<J>1  p 34>  1M  ‘ lt0Po  3p  ^ 1 + ^lR^ 


(28) 


That  is, 
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The  fields  of  all  higher  order  multipoles  will  be  neglected. 

The  boundary  conditions  at  the  surface  of  the  perfectly  conducting 
wire  p = a imply  that, 

AoEzo  * AlEIl  * MlE"l  = 0 

Vjo  * A1E»1  * *¥$!  * « 


(29) 


(30) 

(31) 


on  the  wire.  In  these  equations  Aq,  Aj,  and  Mj  are  unknown  amplitudes  of  the 
zero  order  electric,  first  order  electric,  and  first  order  magnetic  multipoles, 
respectively.  Each  maltipole  field  evaluated  on  the  wire  surface  is  now 
written  in  Fourier  components.  For  example. 


(32) 


Ezo  a Ezoo  + Ezolcos*  * Ezo2cos2^  +*  • 
E”l  = E”llSin*  + E*12sin2*  +*‘- 


Equating  the  coefficients  of  4>  independent,  cosip,  and  sim)>  terms  yield  a 
homogeneous  set  of  three  simultaneous  linear  equations.  Thus,  we  require 
that  the  determinant  of  the  coefficients  be  zero;  that  is. 
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This  represents  the  modal  equation  for  a up  to  first  order  Fourier  components. 
It  remains  to  find  explicit  expressions  for  the  Fourier  components  of  the 
fields,  i.e.,  the  matrix  elements.  Had  we  neglected  the  first  order  terms 

g 

the  equation  would  have  read  Ezoo=0  which  is  equivalent  in  the  thin  wire 
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approximation  to  the  result  of  Wait  and  Chang  and  Olsen 

Representation  of  the  Matrix  Elements 

Recall  that  all  of  the  integrands  involved  here  contain  the  factor 
-UjCkjX+D)  -iXkxy 

e e This  factor  can  be  conveniently  represented  in  the 


following  form. 
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where  <|>o  = -sin  — or  = cos  — . For  our  purposes,  a convenient  form 
of  this  relation  is 
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Upon  substitution  of  this  form  in  the  integrands  it  is  found  expedient 
to  define  two  fundamental  integrals. 
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Now,  all  of  the  integrals  resulting  from  substitution  (35)  may  be  expressed 
in  terms  of  these  two  fundamental  integrals  and  their  derivatives  with 
respect  to  D. 

The  expressions  with  which  we  are  concerned  also  contain  H^^(?kjp') 
which  we  desire  to  write  in  terms  of  p and  <p  ■ This  is  accomplished  by 
the  addition  theorem, 

= l Hm1)(2?D)Jm+n(^klp)ei(m+nH  (38) 
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which  implies  that. 
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Hi(1)UklP')sin  <fr'  = IimHIi1)(2cD)Jm(ck1p)sin  m<j>  (41) 

Substitution  of  (35),  (39),  (40),  and  (41),  and  p * a,  into  (12),  (13), 
(14),  (18),  (27),  and  (29)  yields  the  following  expressions  for  the  matrix 


elements. 
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which  is  certainly  satisfied  when  Ezqo  =0;  as  required  by  Wait  and  by 
2 

Chang  and  Olsen  . It  remains  to  obtain  approximate  expressions  for  the 

fundamental  integrals,  1^  and  I^,  in  order  to  arrive  at  a criterion  for 

validity  of  the  thin  wire  approximation.  We  further  note  that  all  the  I^'s 

and  their  derivatives  vanish  when  the  earth  becomes  perfectly  conducting. 

2 

It  is  then  not  difficult  to  show  that  the  only  acceptable  solution  is  £ = 0 
or  a = ±1.  The  result  then  reduces  to  the  well-known  result  of  a two-wire 
transmission  line. 

Proximity  effect  for  a small  height 

We  now  proceed  to  discuss  the  proximity  effect  as  the  thin- wire 
approaches  the  earth  surface.  In  particular,  we  are  interested  in  the 
case  when  the  height  of  the  wire  is  much  less  than  both  the  free-space 
wavelength  and  the  skin-depth  of  earth  so  that  the  fundamental  integrals 
(36)  and  (37) , and  their  higher  derivatives  can  be  evaluated  in  a manner 
similar  to  that  used  by  Chang  and  Wait^.  Provided  that  |c|d«  1 and 
| n|  d « 1»  the  following  approximate  express  ions  for  lj(a;D)  and  I^^jD) 
hold: 
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where  I and  Ijn  are  independent  of  D and  given  by 
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Successive  differentiation  with  respect  to  D then  gives  the  higher 
derivatives  of  and  1^;  j-2,3,  On  the  other  hand,  the  integration 

on  D on  both  sides  of  (52)  and  (53)  yields  the  result  of 


Va;D)  3 - D{Ho1)^2a,5  - “ [l+lj0(a)]}  (56) 

I2Ca;D)  = I2Q(a)  - D{(-^-)[Ht1)(2CD)  - % - ^ l’Q}  (57) 

1+n 


where  I10(a)  and  I20(a)  are  the  values  Ij  and  I2  evaluated  at 
D = 0.  Thus,  if  now  we  substitute  (52)  - (57)  into  (42)  - (50),  while 
retaining  only  the  leading  terms  in  the  small -argument  expansion  in  D 
and  A,  the  following  approximate  expressions  for  all  the  relevant  field 
components  are  obtained: 


2 2 

b!oo  * T t'2*"  ar  * 

n +1 


where  y » 0.577216 


T f J » 2j » 

30  = 10  _ a 20  ’ 


is  the  Euler's  constant,  A * kjS  and 


2 2 

Ee  ~ — (— ) (r2  i ■> 

Ez01  " it  W ^ 2 * 
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.2  . e A>i  2r_2  , r2Cn 


if  * ($  f?  -i  (-fH} 

n +1 


With  these  approximate  expressions,  the  original  model  equation  as  given 
by  the  determinant  in  (33)  now  takes  up  a much  simplified  form  of 

Etoo  ♦ T sta)  ■ 0 • (< 

or  2 2 

2 2n 

-Cln  * y---  (InCD  ♦ y - i w/2)  ~ i'  ♦ 6(a)  = 0 , « 

A n +1 


where 


;2-i2c2/(n2*l) 


crJ1 


and  the  expression  for  IjQ  = - a I^q  is  given  by  (54)  and  (55). 

Except  for  the  term  6(a),  (67)  and  (68)  now  correspond  exactly  to  the 

approximate  modal  equation  derived  by  Chang  and  Wait5  using  the  so-called 

"thin-wire"  approximation  in  which  the  longitudinal  current  is  assumed 

to  be  angularly  uniform  and  the  azimuthal  current  neglected.  Denoting 

a as  the  solution  of  E®  n(a  ) = 0,  one  then  concludes  that  the  con- 
o zuu  o 

dition  |6(a  )|  « la  3 E®  _(a  )|  must  be  imposed.  We  note  that  a 
1 o o a zuu  o 1 

suitable  solution  of  E200(ao)  = 0 is  known  in  Chang  and  Wait5  to  be 

5o  = ~ l4n  *>l  f*n  2D/A]”1 
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16 


2 2 2 -1 

provided  that  the  condition  |n  |»  zo  » |n  | holds.  Thus,  6(a0)  is 

2 

generally  of  the  order  of  (A/2D)  «,n(2D/A).  We  further  note  that  the 
magnitude  of  first  order  electric  and  magnetic  multipoles  are  related  to 
the  zero-order  electric  multipole  by  (30) -(32).  Using  the  simplified 
expressions  in  (58)  - (66)  and  also  neglecting  terms  of  order  (A/2D)2, 
it  is  not  difficult  to  show  that 


(70) 

(71) 
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large  because  of  the  different  functional  behavior  in  a . In  particular, 

it  has  been  shown  that,  as  a result  of  the  residue  contribution  from  a 

pair  of  poles  in  the  integral  given  in  (37),  both  I3  and  1 4 possess  an 

inverse  square-root  singularity  with  the  branch  points  located  at 
2 i 

<*p  = ±n/(n  +1)  in  the  complex  a-plane.  These  branch  cuts  are  associated 

with  that  part  of  the  continuous  spectrum  which  radiates  along  the  earth 

surface.  Thus,  in  the  neighborhood  of  a * ap,  we  have  from  (37)  and  the 

definition  of  I,  and  I. 

3 4 


T _ 2r  . T 2n2  .2  2,-i  -i2D(n2-l)~* 

l3  s ° V U = TT  («p-«  ) e 
n -1  F 


(73) 


with  the  higher-order  derivatives  and  for  m>0  given  as 


t(»)  _ 

3.4  " 


u(n2-l) 


(74) 


Provided  that  (A/2D)  is  small,  we  can  again  find  the  approximate  field 

expression  from  (42)- (50).  A similar  analysis  then  leads  to  the  con- 

elusion  that  |6(a)|  is  negligible  whenever  |A  yn°|  « | (a  -or)* | . 

2 2 

Thus  except  for  a small  region  near  o»l,  specified  by  |(a  -1)|  « (A/2D y 
according  to  (69),  and  an  even  smaller  region  near  a =ap»  the  thin-wire 
approximation  is  indeed  well  justified. 


Corresponding  results  for  a large  height 

In  the  case  when  the  height  is  greater  than  several  skin-depths  of 
earth  (i.e.,  Inl^2  » 1),  the  fundamental  integrals  (36)  and  (37)  can  be 
approximated  asymptotically  in  a manner  similar  to  that  used  by  Chang  and 
Olsen.  (In  fact  their  W(a)  is  in  their  approximation  just  equal  to  our 
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iirn  I £•)  Beginning  with  1^  we  write 


, r u.+in  -2u,D 

-5 / — e 1 dX 

2 J u. 


itr  (n  -1)  U1 


->r— I [?H(1)(2CD)  + inH^tfCD)] 
n -1  1 ° 


(75) 


where  u2  has  been  approximated  by  -in  because  most  of  the  contribution 
to  the  integral  occurs  where 


r, 2 ^ 2 2 * . r,  X2+a2  , 

u2  = [X  + a -n  ] « -m[l 5—  ] 


2n 


(76) 


X2+a2 


2n 


and  y~ J«l.  As  indicated  by  Chang  and  Olsen,  a better  approximation 


is  u 


. 2 
-in 


2 ~(n2*l)* 


An  approximation  for  I2  can  be  obtained  directly  from  Chang  and 
Olsen  as  mentioned  above.  However,  we  present  here  the  derivation  of  a 
more  approximate  but  simpler  expression.  First,  we  cast  I2  in  the  form. 


. x 
1- 


1 2 a -J-  f H^1)(x)e  nC  e-2iD/n  dx 
n 4 2?D 


(77) 


An  approximation  for  large  2?D  may  now  be  obtained  using  the  large 
argument  asymptotic  expression  for  H^(x).  That  is. 


00 


1 1(1V)X  -2iD/n 

— e e dx 

v6T 


2 CD 


yfl 


n't/TTX- 


erfc 


Itt)/ 
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-2i  D/n 


(78) 


Now,  using  the  large  argument  asymptotic  form  for  the  complementary  error 
function  we  obtain: 


To  the  same  order  in 


1 

1/CD  we  have 


e2UD 

/ccT 


(79) 


2iCD 


yer 


(80) 


To  this  order  of  approximation  g-jy  — *■  2ic  so  that  the  derivatives  of 
the  I's  in  the  matrix  elements  may  all  be  replaced  by  factors  of  2iC  . 

Use  of  (79)  and  (80)  will  result  in  a lowest  order  approximation  for 
large  £0.  if  is  not  very  large  one  must  obtain  more  accurate  ex- 
pressions for  1^  and  as  was  done  by  Chang  and  Olsen.  For  purposes 

of  this  treatment,  however,  it  will  be  assumed  that  CD  is  sufficiently 
large  that  our  asymptotic  forms  are  appropriate. 
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The  approximations  obtained  above  will  now  be  substituted  in  (42)- (50) 
and  the  amplitudes  A1  and  will  be  found  in  terms  of  Aq  and 
parameters.  For  notational  convenience  we  define. 


2UD 


/iir^D 


(81) 


Thus, 


5 -w 


(82) 


i. . .if-msaL ) 

\ 1 + 


2 ~ ' l/(nC); 


(83) 


and 


T • fn-C  (a/n  Q2 
3 Ln-1  1 + 1/  (n£)  J 


l/(nc) 


G = q,G 


(84) 
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n-C  _ (1/n 


>-n  -1  l + l/(n£) 


-] 


g = q4G 


(85) 


Retaining  only  constant  and  linear  terms  in  the  series  for  the  Bessel 
functions  J^Ckjp)  and  evaluating  at  p = a;  i.e.,  k^p  = A,  we  have 
for  large  CD  and  small  £A, 

.2 


Ez00  ~ ^UA)  - (C‘*2iMs)G 


Ezoi  a WK  -2ws>G 


EJoi  = - 2q3)G 
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Using  these  approximations  to  the  matrix  elements  in  (33) , one  then 


obtains  E®0Q  + 6(a)  = 0,  where 


6(a)  = -(^)  (CA)  (C+i2q3)  (C-2qs)  ei4?D  (95) 

The  condition  |6(a)|  is  small  requires  a small  radius-to-height  ratio 


in  addition  to  a small  |ca|.  Letting  be  such  that  *^00 (CQ)  = 


and  defining  A?  to  be  where  is  the  value  of  t,  which 

satisfies  (33)  in  this  approximation  we  find  in  agreement  with  the 
general  theory  of  small  perturbations  that 
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so  that  the  ratio  of  the  first-order  current  component  to  the  uniform 
current  is  given  by 

lJsZl  I i A (ao/n  *o> 


(#)*  1 + 1/  (n^) 


Now,  since  CqD  is  the  large  parameter,  one  might  say  the  currents  are 
perturbed  but  the  propagation  constant  is  not.  The  conclusion  we 
derived  for  a large  height  is  therefore  consistent  with  the  case  of  a 
small  height. 

Conclusion 

The  work  presented  was  primarily  concerned  with  obtaining  a measure 
of  the  significance  of  the  first-order  components  of  the  currents  on  the 
wire  compared  with  the  zero-order  on  the  uniform  current.  Using  a 
multipole  expansion  method,  we  derived  a formal  expression  for  the  modal 
equation  in  the  form  of  a determinant,  i.e.  (33).  Validity  of  the  so-called 
"thin-wire"  approximation  was  discussed  for  two  cases;  one  corresponds 
to  the  case  of  a vanishing  height,  | n | D « 1 and  the  other  for  height 
greater  than  several  wavelengths  of  waves  in  free-space.  In  both  cases 
the  amplitude  of  the  first  order  current  becomes  more  significant  as  the 
height  decreases  before  the  propagation  constant  is  significantly  modified. 
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